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The Mott transition in a two-band Hubbard model involving subbands of different widths is 
studied as a function of temperature using dynamical mean field theory combined with exact 
diagonalization. The phase diagram is shown to exhibit two successive first-order transitions if the 
full Hund's rule coupling is included. In the absence of spin-flip and pair-exchange terms the lower 
transition remains first-order while the upper becomes continuous. 
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The nature of the metal insulator transition in ma- 
terials involving subbands of different widths has been 
intensively debated during the recent years 0, 0, E3 [J IE 
E S H H [H E3- This issue is relevant for the un- 
derstanding of the effect of strong local Coulomb inter- 
actions in systems such as Ca2-a;Sr x Ru04. In this layer 
perovskite the partially filled Ru t 2g bands for x — 2 con- 
sist of wide d xy and narrow d xz ,y Z bands fl3j |. As pointed 
out in 14], the onsite Coulomb energy lies between the 
subband widths: W xz , y z < U < W xy . The usual crite- 
rion with the parameter U /W as a measure of the impor- 
tance of correlations must then be generalized. The pure 
Sr compound is superconducting below 1.5 K 15]. Iso- 
electronic replacement of Sr by Ca leads to an effective 
band narrowing due to octahedral distortions [lf3 ] and a 
metal insulator transition ^7} ■ As a consequence of non- 
cubic crystal fields many other transition metal oxides 
also involve non-equivalent partially occupied subbands. 

A key question in these materials is whether the wide 
and narrow subbands exhibit separate Mott transitions 
or whether single-particle hybridization and inter-orbital 
Coulomb interactions ensure a single transition for all 
bands simultaneously. This issue was studied initially by 
Anisimov et al. |l| and Liebsch 2] for Ca2- a; Sr 2 ;Ru04 us- 
ing simplified band structure models which did not yet 
include the full complexity due to Ca-induced octahe- 
dral distortions. Correlations were treated in the dynam- 
ical mean field theory |l8| (DMFT) combined with the 
non-crossing approximation (NCA) and Quantum Monte 
Carlo (QMC) method, respectively. Since the NCA cal- 
culations (at T = 0) neglected interorbital Coulomb in- 
teractions the results showed separate, 'orbital-selective' 
Mott transitions for the narrow and wide subbands 0. 
In contrast, the QMC calculations (at T = 0.125 eV) in- 
cluded interorbital Coulomb interactions and suggested 
a common transition for all ti g bands Q. 

Recent theoretical studies of the Mott transition in a 
paramagnetic two-band model system have led to appar- 
ent contradictions 0, S S S S H El . Including the 
full Hund's rule coupling Koga et al. |4| found orbital- 
selective metal insulator transitions at T — 0. Neglecting 
spin-flip and pair-exchange terms Liebsch || obtained at 
T > a single first-order transition, followed by a bad- 



metallic/non-Fermi- liquid (NFL) phase. As we argue be- 
low, the various results and conclusions are consistent 
provided that the choice of crucial parameters such as 
temperature and Hund's rule coupling is properly taken 
into account. 

In the present work we combine finite temperature 
DMFT with exact diagonalization |l9j (ED) to deter- 
mine the T/U phase diagram of a two-band system 
consisting of non-hybridizing, half-filled subbands with 
semi-circular density of states of width W\ = 2 eV and 
W% = 4 eV. The subbands interact via intra- and in- 
terorbital Coulomb matrix elements U and U' — U — 2 J, 
where J is the Hund's rule exchange integral. In con- 
trast to the multi-band QMC approach, which includes 
only Ising-likc exchange terms to avoid sign problems at 
low temperatures |20j , ED permits also the consideration 
of spin-flip and pair-exchange interactions. To be specific 
we take J = U/4 which is approximately satisfied in sev- 
eral transition metal oxides. 

The main result of this work is that the T/U phase 
diagram in the presence of the full Hund's rule exchange 
exhibits two successive first- order phase transitions, with 
separate hysteresis loops and coexistence regions. The in- 
termediate region corresponds to the T > analog of the 
orbital-selective Mott (OSM) phase obtained at T = in 
4]. On the other hand, if spin-flip and pair-exchange 
terms are omitted, we find a single first- order transition 
succeeded by a non- Fermi-liquid phase, in agreement with 
previous QMC results 0. Both trends are consistent 
with those obtained by several groups 0, 0, 0, 0] for 
the same two-band model at T = 0. 

The ED/DMFT results are derived from a two-band 
generalization of the approach employed for single-bands 
El I3II3- Since at T > all states of the impurity 
Hamiltonian are used in the construction of the subband 
Green's functions Gi(itu n ), two bath levels per impurity 
level are taken into account (n s = 6 per spin). To check 
the accuracy of this approximation we have evaluated the 
T/U phase diagram of a single band for n s — 3, ... ,6. 
As shown in Fig. 1, the stability boundaries U C \(T) and 
U C 2(T) for n s = 3 are slightly too low. Nevertheless, the 
overall shape of the phase diagram agrees qualitatively 
with the converged results for n s = 6 23]. Thus, we are 
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FIG. 1: Phase diagram for one-band Hubbard model, calcu- 
lated within ED/DMFT for n s = 3, 4, 6. Symbols (x): T = 
transitions obtained in [231. 



confident that in the two-band case n s = 6 also yields 
a reasonable picture of the phase diagram. Preliminary 
results for n s = 8 will also be presented. 

To analyze the metal/insulator transition we study the 
quasiparticle weights Z; = 1/[1 — dRe'Ei(u))/du)\ u= o], 
which in the metallic range can be represented as Zi f=a 
1/[1 — ImEj(io; )/a;o], where Xi(ia;o) is the subband self- 
energy at the first Matsubara frequency. Fig. 2(a) shows 
Zi(U) for J' = J = U /4, where J' denotes spin- flip 
and pair-exchange terms [24|. Two critical regions can 
be identified, each with hysteresis loops characteristic of 
first-order phase transitions. The coexistence areas are 
Uf x {T) < U < U< 2 (T) near 2.0 eV and U> X {T) < U < 
U> 2 {T) near 3.0 eV, where U<(T) and U> n (T) are the 
stability boundaries obtained for increasing (n — 2) and 
decreasing (n = 1) U. Let us denote the true critical 
energies of these transitions as (T) and U> (T) . Be- 
low U<(T) both bands are metallic while above U>(T) 
both are insulating. At the lower transition both bands 
undergo first-order transitions - but in fundamentally 
different ways: Z\{U) becomes very small while Z2(U) 
drops to a finite value. The narrow band therefore under- 
goes a 'complete' metal/insulator transition, whereas the 
wide band exhibits an 'incomplete' transition to a new, 
considerably more correlated phase. This band becomes 
fully insulating near 3.0 eV, where it exhibits a weak sec- 
ond hysteresis loop. To our knowledge this is the first 
time that sequential first-order transitions are identified 
in the T/U phase diagram of a Hubbard model involving 
non-equivalent bands |25j . If such a material could be en- 
countered experimentally, the conductivity as a function 
of pressure would show two consecutive jumps. 

Fig. 2(b) shows Zi(U) for J' = 0, J = C//4, i.e., in the 
absence of spin-flip and pair-exchange terms. The results 
are similar to those in Fig. 2(a), with the important ex- 
ception that the wide band above the lower transition 
is even more correlated and the upper transition is now 
continuous at smaller U |26| . Remarkably, the lower tran- 
sition remains first-order for both subbands. In this case 
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FIG. 2: Zi(U) of nonisotropic two-band Hubbard model, 
calculated within ED/DMFT. (a) J' = J = (7/4, (b) J' = 0, 
J = U/4. Solid (dashed) curves: narrow (wide) band. Results 
for different temperatures are displaced vertically by 0.3. 



the conductivity shows a jump at U<(T) but a change 
of slope at U>(T). The results in Fig. 2(b) confirm the 
picture obtained previously within the QMC for T > 0, 
J' = which showed the existence of a single first-order 
transition followed by a mixed insulating/bad-metallic 
phase [{J |U ■ They also demonstrate that the two- 
band ED/DMFT for n s = 6 is qualitatively accurate [2^| . 

The phase diagrams deduced from the ED results for 
T > 2.5 meV are shown in Fig. 3. For J' = J as well 
as J' — the transition at U<(T) is first-order for both 
subbands. The subsequent transition of the wide band at 
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FIG. 3: Phase diagram for nonisotropic two-band Hubbard 
model, calculated within ED/DMFT. (a) J' = J = [7/4, (b) 
J' = 0, J = (7/4. Solid dots in (a), (b): stability boundaries 
of both subbands near lower first-order transition. Open dots 
in (a): stability boundaries of wide band near upper first- 
order transition. Dashed line in (b): approximate location 
of continuous transition of wide band. Symbols (x): T = 
transitions obtained in |3, 0, Eli (+) : transitions at T = 
31 meV obtained in [H Il2l|. Lines are guides to the eye. 



U>(T) is first-order for J' = J but continuous for J' = 0. 
At U^-{T) the metal/insulator transition is complete only 
for the narrow band. The wide band first undergoes a 
transition to a more strongly correlated phase and be- 
comes truly insulating at the second transition at U^(T). 
The overall shape of the phase diagram for J' = J agrees 
with the one recently obtained by Inaba et al. [2f|. 

Fig. 3 also shows the critical Coulomb energies ob- 
tained at T = 0, i, [13 and T = 31 meV H [Q. As 
illustrated in Fig. 1 for the one-band case, the ED scheme 
with small n s underestimates the critical Coulomb ener- 
gies by about 0.1 to 0.2 eV. Preliminary results for n s = 8 
indicate similar shifts in the two-band case. Taking these 
displacements into account the ED results shown in Fig. 3 
are in excellent correspondence with those obtained at 
T = and T = 31 meV. Evidently the conflicting con- 
clusions reached in [mSIM&IEljIsOA Ui| concerning 
the nature of the Mott transition in multi-band systems 
were caused by different behaviors obtained for T = 
vs. T > and J' = J vs. J' = 0. Accounting for these 
different parameter choices, the DMFT treatments are 
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FIG. 4: Self-energy of wide band for J' — J in intermediate 
phase at U = 2.4 eV. (a) n s = 6; (b) n 3 = 8. Solid dots: 
T — 10 meV; empty dots: T = 5 meV; (+): T = 2.5 meV. 



consistent. 

Quasiparticle spectra derived within the QMC /DMFT 
for J' = showed that in the intermediate phase the 
self-energy of the wide band at small uj n deviates signifi- 
cantly from metallic ~ lu h behavior |30j . Accordingly, the 
quasi-particle spectra show a pseudogap which for larger 
U gets more pronounced |3ll | , indicating progressive non- 
Fermi-liquid properties. A true gap opens at U^(T). 
Thus, Z 2 (U) > in the region U<(T) < U < U>(T) does 
not imply existence of quasiparticles. As indicated in 
Fig. 3(b) the intermediate 'orbital-selective Mott' phase 
in the absence of spin-flip and pair-exchange is in fact a 
mixed insulating/NFL phase. The same trend is found 
using the T > ED/DMFT 32]. 

For J' — J, the low-frequency analysis of ^(iuin) is 
more intricate. As shown in Fig. 4(a) for n s = 6, T,2(iw n ) 
reveals deviations from metallic ~ oj n variation, giving 
rise to small pseudogaps in the quasiparticle spectra. 
This behavior is incompatible with ImE2(w) ~ u> 2 and 
suggests that, as for J' = 0, in the intermediate phase the 
wide band at T > violates Fermi-liquid theory. The ex- 
tension of the present ED approach to n s — 8 indicates, 
however, that the additional bath levels are important for 
the low- frequency variation of T,2(iw n )- As can be seen 
in Fig. 4(b), the deviations are absent and the ~ uj n vari- 
ation is consistent with Fermi-liquid behavior. In fact, 
the shoulder near u>o = 0.06 suggests that Fermi-liquid 
properties persist up to about T = ujq/tt w 20 meV. 
Thus, the OSM phase in Fig. 3(a) is the T > analog of 
the orbital-selective Mott phase identified first by Koga 
el al. at T = 0. A more complete discussion of the 
results for n s — 8 will be given elsewhere • Because of 
finite size limitations of the present ED/DMFT scheme 
associated with the small number of bath levels, a pre- 
cise determination of low-temperature properties is not 
possible. Nevertheless, an approximate extrapolation of 
£2(io> n ) indicates that the T — > limit for J' = J sat- 
isfies Fermi-liquid criteria, in agreement with previous 
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In view of the importance of spin-flip and pair- 
exchange terms for the Mott transition in multi-band 
materials [33j it is desirable to investigate the T/U phase 
diagram by methods which permit adequate treatment of 
the complete Hund's rule matrix, for instance, a two- 
band extension of T > numerical renormalization 
group studies [^i see also |2^|. Because of sign prob- 
lems, recent QMC extensions including spin-flip andpair- 
exchange terms are limited to T > 1/6 eV 3> T c |]| or 
T = |loj|. Also, other recent works [1,0 employing a 
variety of quantum impurity methods deal so far mainly 
with T = and do not yet allow the identification of the 
multi-band Mott transition at general T, U values. 

In summary, the T/U phase diagram of the Hubbard 
model involving half-filled, non-equivalent subbands is 
shown to be remarkably rich. The competing kinetic 
energy scales, coupled via Coulomb and exchange ener- 
gies, give rise to sequential first-order phase transitions. 
The lower transition separates a purely metallic phase 
from a mixed phase with insulating narrow and strongly- 
correlated wide subbands. The wide band becomes in- 
sulating at the second first-order transition. Omission of 
spin-flip and pair-exchange terms enhances the correla- 
tions in the wide band in the intermediate phase so that it 
no longer satisfies Fermi-liquid criteria, and modifies the 
upper phase transition from first-order to continuous. 

For the analysis of experimental data of materials such 
as Ca2- a; Sr 2: Ru04 it is necessary to account also for hy- 
bridization between orbitals. Preliminary studies of this 
effect within two-band models for T ^S> T c Q and T = 
pj.lll| su gges t significant changes. Moreover, spatial fluc- 
tuations |35j and deviations from half-filling might play 
a decisive role close to the Mott transition. More work 
is needed to investigate whether both first-order transi- 
tions persist in the presence of these effects, or whether 
the weak first-order behavior of the upper transition dis- 
appears and only the dominant lower transition survives 
as the common first-order Mott transition for all bands. 
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